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1 Introduction 

Our approach to problems of rarefied gas dynamics is based on direct so¬ 
lution of the complete kinetic Boltzmann equation. The main difficulty of 
solving the Boltzmann equation is the calculation of the multi-dimensional 
collision integral. In this paper, we apply our generalization of the conser¬ 
vative discrete ordinate method of Tcheremissine |1 (which was originally 
developed for a single gas) for binary gas mixtures in the case of cylindrical 
symmetry (Raines A. A.) |2|. For the evaluation of collision integrals we 
use the projection method which ensures the strict conservation of mass, 
momentum and energy. The conservativeness of the method is achieved by 
a special projection of the integrand values, calculated at non-node points, 
to nodes of the momentum grid that are closest to them. 

Using this method, we now solve the problem of interaction of a two- 
component supersonic jet with a normally posed wall. We discover the effects 
of the inflow of gas on a cool wall with the mirror and diffuse reflection laws 
and the appearance of the wall Knudsen layer. The results are compared 
with the gas-dynamical solution and with the papers dealing with a similar 
problem for a single gas |3,[4j. We obtain a good agreement with all those 
results. 

2 Description of the method 

The system of the Boltzmann equations in the momentum space for two 
gas components consisting of hard sphere molecules with the masses m, and 
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diameters di reads 
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where f t is the distribution function which depends on the vector of momen¬ 
tum p, the vector of configuration space x and the time t. 

Direct and reverse collision integrals have the form 
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In the kinetic momentum space we have the following relations between the 
vectors of momentum before and after the collision 
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n = n(cos(9, sin 0 cosy?, sin 0 sin 99 ). 

Here n is a unit vector directed along the interaction line of molecules, g Jt 
is their relative velocity, 6 and p are collision angles. 

We introduce in the limited domain of the configuration space a fixed 
grid. We impose limits on the momentum variables in ([T]) - ([3]) by introduc¬ 
ing a domain Q of volume V. In H we construct a discrete grid containing 
Nq equidistant points pp with the step h which results in the discretization 
of the distribution functions and collision integrals as follows 
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The Boltzmann equation in a discrete form becomes 
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dt mi dx 

Equation Q is solved by the splitting procedure. On each interval At 
we split the process into the two stages, free-molecular flow and collisional 
relaxation described by the following equations 
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Let us consider the integral operator 
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Taking for <J> a three-dimensional 5-function, we reduce the collison integrals 
to the form 

L i(p*) = ~P) + Kp* -P*)) ( 9 ) 
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Transforming the integral operator to cylindrical coordinates and taking into 
account the relations p = p(p, p, 7 ), dp = pdpdpd'y we obtain 
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Introduce the uniform integration grid: p au , p av , pp u , pp u , 7 „, 7 *^, 9 U , p v 
with N v nodes. The multiple integral is calculated as the 8 -fold sum over 
all the nodes while the distribution functions do not depend on 7 ^, 7 *„: 
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Combining ([6]) and ( |11[ ) we obtain 
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We replace the expressions in parentheses in (|12j) by using the relations 
${P* - PLJ + KP* - Pf)J = r ?i 6 (P* - PK+s) + S(P* - P^+g)) ( u ) 


The coefficients r? can be found from the conditions of conservation of the 


density, kinetic momentum and energy in the decomposition (14) for a pair 
of cells including their vertices. For economy of computations it would be 
preferable to have a decomposition with a minimum number of terms, so 
that expression (14) becomes 
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If we use © in © and combine with 0 , then we obtain 
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The expressions (13) and (16) define the conservative discrete ordinate 
method if the coefficients have been already found. 


3 Formulation of the problem 

Consider the problem of the inflow of a binary gas mixture upon a wall with 
mirror and diffusive laws of reflection of the gas from the wall. The half-space 
x > 0 is being filled with the gas moving with velocity U (or momentum 
P ) in the negative direction of the x-axis. The density and temperature of 
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the gas are equal to no and To, respectively. At the instant t = 0 we set an 
immovable wall at the point x = 0. At t > 0 the interaction of the gas with 
the wall starts. The system of Boltzmann equations in momentum space for 
this problem for the two gas components consisting of hard sphere molecules 
with masses rrii and diameters di has the form 0- 

A solution of the problem depends on four variables, so that the distri¬ 
bution functions have the form = fi(t,x,p, p). As an initial condition for 
the distribution functions we take Maxwell functions with parameters of the 
unperturbed flow 
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where n^o is the number density of ith component of the mixture, no = 
nio + n 2 o, m is mass of the mixture. For the calculations, an infinite domain 
of the variation of x is replaced by the segment [0, L\ and at x = L we 
set up conditions © for molecules flying into the domain, which yields a 
boundary condition on the free boundary surface 
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On the left boundary of domain [0, L\ we set up conditions for p ^ 0. In the 
case of mirror reflection from the wall, the boundary conditions have the 
form 

limit = 0 ,x,p ^ 0 ,p) = fi(t,X , -p,p). (19) 

For diffusive reflection with the full accommodation the distribution func¬ 
tions of reflected molecules are assumed to be maxwellian with the temper¬ 
ature equal to the wall temperature: 
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Here T w is the wall temperature, rii W is the density of reflected molecules 
which is found from the condition of non-percolation: 

p< o 
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Thus, the problem of the inflow of a rarefied gas on a wall is posed for 
system 0 with initial conditions boundary conditions (|18|) on the 


free boundary surface and boundary conditions (19), 


on a solid wall. 


Graphical representation of solution results can be seen in figures [l}|6j 


4 Dimensionless quantities 

Here we introduce the dimensionless quantities corresponding to all physical 
variables which are involved in our solution method. 
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5 Parameters of calculations 

All calculations were performed at the laptop computer Sony VAIO, proces¬ 
sor Intel(R) Core(TM) 2CPU, 1.66GHz + 1.66GHz, 1.00GB of RAM. 

Time for an iteration step is 2.4 sec., time for the whole problem until 
t = 30t is 120min (r is the mean free time). 

Numerical parameters of calculations: 

M = 2.8773, U = -2, m 2 /m 1 = 0.5, 

n 2 Jn lo = 0.5/0.5, d 2 / d\ = 1, Kn = 1. 

The grids used for calculations: 

Npp = 1250 = (25 + 25)25 - the number of nodes of the momentum grid 
with the step h = 0.2 

NX = 180 - the number of nodes of the x-grid with the step h x = 0.2. 
NOD = 66000 integration nodes, DT = 0.01 is the time step. 
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6 Comparison of numerical results with the gas- 
dynamical solution 

According to gas dynamics equations, a domain is created with constant 
values of macroscopic parameters behind the reflected shock wave. For pa¬ 
rameters specified in section 5 this gas-dynamical solution is the following 

rigd = 2.936, nigd = n 2gd = 1.468, T gd = 3.4396 

Here n g d , n \ gd , n 2gd are numerical densities of the mixture, first component, 
and second component, respectively, while T gd is the temperature of the 
mixture. 

Our results obtained by the direct numerical solution of the Boltzmann 
equation are described below. 

1. Mirror reflection. 

When T wa u = T gas = 1, there is the zone of constant values for density and 
temperature 


n = 2.93598, m = 1.48607, n 2 = 1.451 
T = 3.44089, T\ = 3.44135, T 2 = 3.4375. 

2. Diffusive reflection. 

a) When T wa u = 0.5, T gas = 1 and t > 18r, x > 12.5A, we have the region 
of constant values for density and temperature 

n = 2.9348, T = 3.3435. 

b) When T wa u = T gas = 1 and t > 18r, x > 12.5A we have the region of 
constant values for density and temperature 

n = 2.9313, T = 3.3135. 

We observe a good agreement of our results with the gas-dynamical solution. 

7 Comparison of our results with the results ob¬ 
tained by other methods 

Among other publications studying this problem, we should mention two 
earlier works that restrict themselves to a one-component gas |3,4j. In [3|, a 
conservative difference scheme is applied for solving the Boltzmann kinetic 
equation by the splitting method, with a good agreement of results shown 
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in Figure [Tj In |]4i] a scheme of the second order of accuracy with respect 
to At is applied on the basis of the model Krook equation. For the sake 
of comparison with |4| we have solved the problem of reflection of a shock 
wave from a wall with a reasonable agreement of our result with [4] shown 
in Figure [8j Initial test of the projection method was made on the problem 
of a shock wave in a binary gas mixture where a comparison of results was 
made with the work |5| for a wide range of parameters. An example of such 
comparison is shown in Figure [9j 
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Figure 1: Mirror reflection. Profiles of temperatures of the mixture compo¬ 
nents at different times. 
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Figure 2: Diffuse reflection. Profiles of numerical densities of the mixture 
components at different times. Temperature of the wall is equal to 0.5. 
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Figure 3: Diffuse reflection. Profiles of numerical densities of the mixture 
components at different times. Temperature of the wall is equal to 1. 
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Figure 4: Diffuse reflection. Profiles of velocities of the mixture components 
at different times. Temperature of the wall is equal to 1. 
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Figure 5: Diffuse reflection. Profiles of temperatures of the mixture compo¬ 
nents at different times. Temperature of the wall is equal to 1. 
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Figure 6: Diffuse reflection. Profiles of the heat flow of the gas mixture at 
different times. Temperature of the wall is equal to 0.5. 
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Figure 7: Comparison of results of Raines (Raynes) with the paper |3j for 
temperature profiles for a one-component gas. 
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Figure 8: The problem of reflection of a shock wave from a wall. Comparison 
of results of Raines (Raynes) with the paper |4j for density profiles for a one- 
component gas. 
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Figure 9: The problem of a shock wave in a binary gas mixture. Comparison 
with results of j5|. Profiles of density, velocity and temperature for the 
mixture. 
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